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Abstract. We present a new six-parameter family of potentials whose solutions are 
expressed in terms of the hypergeometric functions 3F2, 2F2 and 1F2. Both the 
scattering data and the bound states of these potentials are explicitly computed 
and the peculiar properties of the discrete spectrum are depicted in a suitable 
phase diagram. Our starting point is a third-order formal eigenvalue equation of 
the hypergeometric type (with a suitable solution known) which is transformed to 
, the Schrodinger equation by applying the reduction of order technique as the crucial 

s ! ' first step. The general preconditions allowing for the reduction to Schrodinger form 

of an arbitrary eigenvalue equation of higher order, are discussed at the end of the 
article, and two universal features of the potentials arising this way are also stated 
and discussed. In this general scheme the Natanzon potentials are the simplest 
special case, those presented here the next ones, and so on for potentials arising 
C^") ■ from equations of fourth or higher order. 
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Qlh| I. INTRODUCTION 

The search for exactly solvable quantum mechanical potentials has a long history 
dating back to the early days of quantum mechanics. In 1971, Natanzon 1 appeared 
to have closed the subject by explicitly constructing the full set of potentials V(x) 
whose solutions can be expressed in terms of the hypergeometric functions 2-F1 and 

The subject received a new impetus by the addition to the available tools for 
exact solution of the Darboux transformation 2 : A mapping of the Schrodinger equa- 
tion to itself which can be applied to any solvable potential and produce an infinite 
\& \ chain of new ones, usually with predetermined spectral changes with respect to the 

initial potential (addition or removal of bound states, etc.). 
\Q [ Given though that the solutions of these "derived potentials" are also expressed 

in terms of the hypergeometric functions 2F1 and 1-F1 - or, to be precise, in terms 
of linear combinations of these functions and their derivatives 3 - it is reasonable 
to view these potentials as simple "derivatives" of the initial ones. So, a complete 
description of the now-known exactly solvable potentials is, we believe, this: Exactly 
solvable are the Natanzon potentials and their Darboux derivatives. 

In view of the above, it is clear that new solvable potentials can be discovered 
only if we enlarge the set of functions within which their solution is sought for. 
And the most natural such enlargement is the full set of hypergeometric functions 
pF q . These are certainly equally exact functions as %Fi or 1F1, the only difference 
being that they satisfy linear differential equations of order higher than second. 
Specifically, the function p F q satisfies an equation of order q + 1 where for any 
given q the index p takes the values ^ p ^ q + 1 characterizing the various types 
of hypergeometric functions of the given order. 

In particular, if we restrict ourselves to hypergeometric functions of the third 
order -e.g., the functions 3F2, 2-F2, 1-F2 and 0F2- then it is clear that their "use" for 
constructing new exactly solvable potentials presupposes, as a necessary first step, 
the reduction of order of the third-order equation satisfied by these functions. But it 
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is also necessary that this initial third-order equation has an eigenvalue parameter 
A suitably located so that it takes the position of energy in the final Schrodinger 
equation produced by the reduction. How this can be done will be examined in the 
following section. 



II. THIRD-ORDER EIGENVALUE EQUATIONS REDUCIBLE TO THE SCHRODINGER 
FORM. 

As is well-known, if a particular solution y = yo of a linear differential equation 
is available, the transformation y — yoY -that is, the factoring out of the known 
solution- eliminates the Y term of the new equation and so the further substitution 
Y' = u lowers its order by one. In the case of a linear equation of order three 

a{x)y"' + b(x)y" + c{x)y' + d{x)y = (1) 
the application of this procedure leads to the second order equation 

u >> + (l +3 yk) u > + (i +2 b -yA+3 y A) u =o, ( 2 ) 

\a yo) \a a y y J 
where u = (y/yo)' ■ With the further substitution 

u ( x ) = Vo 3/2 ex P (~\ J ~^ dx ^J U ^ ( 3 ) 
equation (f2]) is transformed into the canonical form 



ya 2\aJ 4\aJ 2 a y 2 \y a J 4 \y J J 

where the first derivative term is missing. From (fj]) it is now clear that it can be 
cast into the so-called Liouville form 

U" + (Xw(x) - v{x))U = (5) 

-which is readily reduced to the Schrodinger equation- if the following conditions 
are met: i) The coefficients a(x), b(x) and the special solution y = yo of (JTJ) , are 
independent of the eigenvalue parameter A, ii) the coefficient c depends linearly on 
A. That is, 



c(x, A) = cq(x) + Aci (x). (6) 

If these conditions are met then Q is indeed a Liouville equation with weight 
function w(x) — c\{x)/a{x) and Liouville potential -this is a suitable name for the 
function v(x) in ([5])- as follows 

v[x) = _££ 1 fbV 1 (b_V _ 1 ^ _ 8 (yjA' _ 8 (yj,y 
V ^ a 2\aJ 4\aJ 2 a y 2 \y J 4 \y J V ; 
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As far as equation ([5]) is concerned it is well known that it can be reduced to the 
Schrddinger form 

^"(z) + (\-V(z))^(z)=0 (8) 

-where primes stand now for derivatives with respect to z— with the so-called 
Liouville transformation 

i,{z) = i>{x{z))=w 1 ' i U{x)\ x=x{z) (9) 
where the final position variable z is defined by 



z'(x) = \/w(x) => z(x) = / \Jw{x) dx (10) 



and the potential V(x) -from now on we will always assume that x — x(z)- is given 

by 

= „(») + {z,s} 
w(x) 

where {z,x} is the so-called Schwartzian derivative -or simply the Schwartzian- of 
the function z(x) defined by 

i / II \ 1 i / II \ 2 i / I \ ' 1 / / \ 2 

1 / z \ 1 / z \ 1 I w \ 1 / w 



<*•*> = 2 -lUJ =iUJ -leUJ • (12) 

It is now clear from (fT2"]) that the integration constant in (fluT) is eliminated and can 
be ignored. Concerning the requirement that the special solution yo is independent 
of A -in spite of the fact that A appears in the equation- this leads also to the 
condition that the coefficient d in ([T|) is linearly dependent on A like c. That is, 

d(x,X) = d (x) + Xd 1 (x), (13) 

whence equation ([TJ can finally be written as 

(L + XM)y = 0, (14) 

where the operator L has the full order of the equation -that is, three - while M is 
of first order. We claim therefore that any formal -i.e., with no boundary conditions 
specified- eigenvalue equation of the form (|14p can be reduced to the Schrodinger 
equation, provided that the solution y ~ ya(x) of the first-order equation My = 
is also a solution of Ly = 0, and therefore also of the complete equation (fT4"|) 
for all values of the parameter A. A suitable generalization of this proposition for 
equations of arbitrary order will be given in section VII. 



III. EXACTLY SOLVABLE EIGENVALUE EQUATIONS OF THIRD ORDER AND THE AS- 
SOCIATED POTENTIALS 

The next step is obvious. Starting with the most general third-order linear equa- 
tion, solvable in terms of hypergeometric functions, we write it in the form (| 14[) 
-with its parameters restricted so that the solution of My = is also a solution 
of Ly = 0- and then transform it into Schrodinger form producing along the way 
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the respective solvable potential. The most general third-order equation we need is 
written as 



x 2 



(1 + ux)y"' + x(a + j3x)y" + (7 + Sx)y' + (~ + c) y = (15) 



whose form is readily recognizable if we introduce the concept of dimension d = 
m — n of the typical term x m y ljL ^ (of an arbitrary linear equation with polynomial 
coefficients) and say that (|15[) is a bidimensional equation. That is, an equation 
whose terms can be grouped in only two different sets, each set having a definite 
dimension. Bidimensional equations are important because their solution in power 
series leads to a two-term recursion formula -e.g., ak+t = /(fc)afc- that allows us 
to determine the general series coefficient in closed form as a function of k. 

Thus the series represents an exactly known function and the pertinent equation 
can be classified as exactly solvable on purely mathematical grounds. Furthermore, 
since the dimension d of the typical term x n y^ n ' is but the displacement on the 
exponent of a typical power x k brought about by the operator L — x m d n , it 
follows that for a bidimensional equation with dimensions d\ and d 2 , the integer £ 
in the recursion relation a^+t — /(fc)a& -which is clearly equal to the difference of 
the two displacements experienced by the general power of the series-solution when 
substituted into the equation- will be I — d 2 — d\ (d 2 > d\ by convention). The 
quantity I will be called the step of the equation since it indeed tells us that the 
series-solution proceeds in steps of size £. Based on the preceding discussion, Eq. 
(| 15(1 is uniquely specified by the statement that: It is the most general bidimensional 
equation of third order and step unity. 

Concerning our assertion that Eq. p"5)) is solvable in terms of hypergeometric 
functions, this follows directly from the existence of a two-term recursion relation 
and will be discussed further in Appendix A where we will also provide a simple 
recipe for arriving at the solution using only elementary algebraic operations. Note 
also that the exponents of power behaviors at zero and infinity, for any bidimen- 
sional equation, are determined by the conditions 

L lX » = L 1 {n)x>* +dl = 0, L 2 x u = L 2 (v)x v+d * = (16) 

i.e., from the roots of the characteristic polynomials Li(fj.) and L 2 (y) of the uni- 
dimensional operators L\ and L 2 representing the components with the least and 
the largest- dimension respectively, of the bidimensional operator L = L\ + L 2 . 

From the preceding discussion it follows that our starting point should be the 
solvable, third-order eigenvalue equation 

x 2 (l + ujx)y"' + x(a + (3x)y" + (7 + 8x)y' + (- + c) y+ 

\({f+5x)y' + U+CU) =0 (17) 
which has the same form as eq. with 

L = x 2 (l + ojx)d 3 + x(a + (3x)d 2 + (7 + 5x)d + ( J + C) , (18) 



and 



x 



M = (7 + 8x)d + ( - + C ) {d = d/dx). (19) 
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where 7, S, e and £ are new, arbitrary parameters, independent of 7, 5, e and 
£ of Eq. (fT8|) but with the same role as these, hence the similar notation. With 
M given by the expression above, the general solution of My — takes the form 
V = Ho — x ' s ( x + p) s -with suitably chosen s, s' and p— but here we will restrict our 
discussion to the case s' = 1 since then y$ can certainly be also a solution of Ly = 
as it is a terminating series about x — and the equations of hypergeometric type 
are especially these admitting such solutions. 

By regarding now s and p as given parameters and requesting that yo — x s (x+p) 
be a simultaneous solution of the equations My = and Ly = 0, we find that the 
operator M can only have the form 

while L has the form of (|18j) . with the parameters 5, e and £ no longer independent, 
but expressed in terms of a, /?, w, p and s as follows: 

^ = — (3(1 — ujp)s 2 + (2a - 2f3p + 3ujp - 3)s) 

P 

e = -(s 3 + {a- 3)s 2 + (2 - a) a) (21) 

C = -- ((3 - 2up)s 3 + (2a - f3p)s 2 + (2a ~(3p + 2up - 3)s) . 
p y 

As for the parameter 7 it can be set equal to zero since it only appears as an 
additive constant to the potential and can therefore be ignored. 

Our assertion becomes now quite specific: We claim that the formal eigenvalue 
equation of third order 

x 2 (l+ux)y'"+x{a+Px)y"+ (sx + A (l + ~Y\ y'+ — + ( - V = 

9 . (22) 

-with (5, £ and C as in ()2ip - is reduced to a Schrodinger equation which must 
therefore be solvable, since Eq. (|22|) is solvable. More specifically, the expression 
w(x) — ci(x) / a(x) implies that the weight function in the relevant Liouville equa- 
tion (Eq. ([5])) will be equal to 

, x ^° In x =4> x = e z 

V ' \x^oo 1 lna , ^ x = &s /pUjz 

where we have assumed p and u> to be positive so that the change of variables 
z = z(x) maps the region < x < 00 -between the singular points x — and 
x = 00 of the initial equation- into the full region — 00 < z < 00 of the final 
position variable z. As for the complete expression of the function z = z(x), this is 
a simple elementary function which, neverthelles, does not seem to be invertible in 
terms of elementary functions even though the existence of an inverse is guaranteed 
by the fact that z'(x) > in the region of positive x. But just as in the case of 
Natanzon potentials, this poses no problem in solving the Schrodinger equation 
and calculating its spectral properties, since only the asymptotic forms of x(z) are 
involved in these calculations and these are explicitly known (Eq. (|23[) ). We note 
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also that the initial variable x is much more suitable for working out the solution; 
so we will use this variable henceforth in the understanding that at the end one 
has to make the substitution x x{z). 

Based on the preceding discussion, we can now readily construct the potentials 
derived from (j2"2"|) by simply applying formulas (UJ through (|12p . Let us begin with 
the Schwartzian term V s (x) = {z, x}/w(x) which is given by 

_ p 4w 2 x 4 + (4w + \2u?p)x z + (3 + 18wp + 3u 2 p 2 )x 2 + {Up + Aujp 2 )x + Ap 2 
si - X > ~ 16 (x + p) 3 (l + wx) 

(24) 

while for the full potential V(x) — (v(x)/w(x)) + V s (x) we get 
... , Ax 4 + Bx 3 + Cx 2 +Dx + E 

V(X) = -, r^— r (25) 

V ' (x + p) 3 (l+ux) y ' 



whereby 



9 „,. oY..3 , {Zo„ o_. o... , oY.. , (PP 



A = ( -poj — 3 j ujs + y^PP ^ 2a — 3u>p + 3 j ws + ^— h pu> — j3p ) ui 

B= ('^V-^P-3) S 2 + 

— —oj 2 p 2 + 3uip + 3uip 2 (3 + —p(3 — —Oicop — 2a 
— (p 2 (3 2 — 3ojp 2 f3 + pa/3 — 3aup + 3uip — p/3) + -ujp(3ujp + 1) 

C = | (^9{uj 2 p 2 - 3)s 2 + 6(4 - 3a - 2aujp + 2(3p + ujp 2 (3 - u; 2 p 2 )s (26) 

+p 2 /3 2 - 2ujp 2 l3 + Aa/3p - 2/3p ~ Watop + a 2 - 4a 
3 2 2 9 15 

+ -L0 P +^P+~ 
D 2 

D = ^ (6(cjp - 3)s 2 + (6/3/0 - 2acj/9 - 12a + 18) s 

+(2a/3p - Aauip + 2a 2 - 6a + up + 3)) 
E= ^|s 2 -i(a-3) S +i(a-l) 2 )p 3 . 

Note that even though the potential is a simple rational junction of the initial 
variable x, it nevertheless depends in a complicated way on the five dimensionless 
parameters a, /3, to, p and s of the problem at hand. (There are five and not six 
parameters as we had initially claimed, since the sixth is a scaling parameter of 
z which we implicitly set equal to unity -see, e.g., the first numerical coefficient 
in (|2"0|) - so that, together with the substitutions h = 2m = 1, we can arrive at a 
complete system of units for the Schrodinger equation.) The form of V{z) -which 
is but the form of V(x) in the interval < x < oo stretched to — oo < z < +oo- 
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will typically be as shown in Figure 1 with Vb (= E/p 3 ) and Voo (= A/ui) given by 

/9 \ /3 \ /3 2 p ^ 
Voc — ( jpw — 3 J s 2 + ( -^p — 2a — Sup + 3 j s + — h pw — /3p. 

Even though the expression for the potential in terms of its parameters is quite 
complex, the solutions of the respective Schrodinger equation look much simpler! 
To a large extent, this simplification is due to a different parametrization of the 
problem which occurs naturally once we realize that Eq. (f22j) has solution 
provided A takes the value 

A = -(3.s 2 + (2a- 3)s) (28) 




so that the corresponding eigenfunction ipo(x) is 



^»(») = — (^F 3 — ' (29) 



where q and r are given by the expressions 



3s + a- l 3.s + Qg/cj - 2 
9= ~ 2 , r = . (30) 

But since for x — > and i->cowe have ipo{%) ~^ % q and ipo{x) — > x r respectively, 
it follows that (|2U)) will satisfy the boundary conditions at the origin and at infinity 
-i.e., at — oo and +oo of the variable z— only if q > and r < 0. In this case, 
ipo{x) represents the system's ground state since it has no nodes in this region. The 
crucial finding now is this: If we take the value (f2"5|) as the reference level for the 
potential energies and eigenvalues of the problem, using instead of V and A the 
quantities 

U = V-\ Q , e = A-A (31) 

then a remarkable simplification ensues. U and e do not depend separately on the 
parameters s, a and /3, but only on their combinations q and r. In other words, when 
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we make this (parametrically dependent) change of reference level, the number of 
(dimensionless) parameters of the problem gets reduced from five to four ; q, r, 
p and oj. Note, for example, that subtracting Ao from the asymptotic values (|27[) 
gives 

U = V -X = q 2 , Uoo = Voo - Ao = pur 2 , (32) 

which are much simpler expressions than those before, with the parameters now 
being q,r,p and oj. Let's call this way of parametrizing the problem invariant 
parametrization; we shall work with this from now on. The expression for the po- 
tential U(x) now becomes as in (|25[) . i.e., 

TTI , Ax A + Bx 3 + Cx 2 + Vx + £ 

U{x) = r (33) 

(x + py(i + lox) 

but with new coefficients A, B, C etc., which are given by 
A = u 2 pr 2 

B = 2uj 2 p 2 r 2 + 2ujpqr - iopq + uj 2 p 2 r + ^wp(l - top) 

C= p(q 2 + uj 2 p 2 r 2 + Aujpqr - (1 + ojp)q + up(l + ojp)r + A(l - wp ) 2 ^ (34) 

V = p 2 ^2q 2 + 2uipqr — q + uipr + ~^(y>p — 1) 
£ = p 3 q 2 . 



IV. RESULTS: THE FULL HYPERGEOMETRIC CASE 

As explained in Appendix A, if we use the above parametrization, the solution 
tp(x) of the Schrodinger equation satisfying the boundary condition ip(x = 0) = 
ip(z = — 00) =0 in the discrete spectrum region -or tp(z — > —00) ~ e lkz in the case 
of continuous spectrum- can be written as 



^ {X) = (x + p)V* \ X{X + P)F ' + + P(a + (35) 

where F =3.^2 (a, b, c;d,e; —lux) and a, b, c, d, e are given by 



Q + vV - e d=2-q+^/q 2 



b = ^q 2 -e+ ^r 2 -ge-a e=l+2yY 



(36) 



VT-e - \/r 2 - ge-a 



and 



a = q-r-l, g=—- (37) 

pLJ 
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In fact it can be shown that expression (|35j) can be also written in the equivalent 
and more explicit form 



= j^TpF 1 x 

x (p(a + 1) 3F 2 (a, b, c; a + 1, e; — uix) + ax s F 2 (a + 1, b, c; a + 2, e; — oox)) (38) 

where no derivatives of 3.F2 enter. To the best of my knowledge expressions like (|38[) 
-with higher-order hypergeometric functions present in the solution of a Schrodinger 
equation- appear for the first time in the literature. 

To impose the boundary condition at x — > +00 (i.e., at z — > +00) we also need 
the asymptotic relation 

r(d)r(e) r(fr-a)r( c -q) r(#) r(o-6)r( c -6) 6 

* [X> ™ r(6)r(c) T(d - a)T{e - a) [ ' + T(a)T(c) T(d - b)T{e - b) 1 ' 

T{d)T{e) T{a-c)T{b-c) 

+ f(o)f(6) r(d-c)r(e-c) {UJX) (39) 

in conjuction with the fact that the first-order operator M. acting on F 
-Eq. ([33]) - annihilates the term x~ a , for the same reason that M annihilates also 
the solution yo- Thus, we come to the following conclusions. 

A. Scattering states 

Denoting by r#(fc, k') and r^{k, k') the reflection amplitudes from the right and 
left respectively, we have 

r + i^gk' Y{2i^gk') 



r R (k,k') 



- - i^/gk' T(-2z^gk') 

T(-a - ijk + ^gk'))T(a + 1 - »(fc + ^gk')) ^ 
r(-a - i(k - ^/gk'))T(a + 1 - i(k - y/gk')) 



(40) 



-q + zfc r(i + 2 ? fc) 

rL (k,k) = x — — x 

q + ik I (1 — 2zk) 

/ r(-o- - ijk + ^gk'))T(a + 1 - i{k + VfffcQ) _ 2ifc 

rX-a + l (fc - vp^rv + 1 + ^-^')) ^ ^ J 

where fc and k' are the wavenumbers for z — > —00 and z — > +00 respectively. 
In other words, we have k = ye — q 2 , k' = ye — (r 2 /g). For the corresponding 
reflection probability P r -which is independent of the direction of incidence- we 
get 

cosh 2 ir(k — yfqk') — cos 2 7T(T , % 

P P (fc, fc') = • ( 42 ) 

cosh 7r(fc + yfgk') — cos 2 7rcr 

B. Bound states 

In this case, the requirement for vanishing ip( x ) a t infinity leads -due to (|39|) - 
to the three conditions i) a — —n, ii) b = —n and iii) e — c = —n. From these we 
obtain (see Appendix A for details) the respective regions, red, green and blue, of 
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the phase diagram in Figure 2. More specifically, the condition a = —n is actually 
restricted to the value n — and has a respective (unique) eigenvalue e = 0; while 
b = —n and e — c = — n lead to the equations (I) and (II) for the green and blue 
region respectively. 

Note however that in order for solutions of, say, Eq. (I) to exist, it is not sufficient 
that a > 0, but we must also have a > f m i n , where / m ; n is the minimum value 
of the function /(e) = ^/ q 2 — e + \/ r 2 — ge, in the left hand side of (I), which is 
monotonically decreasing and will therefore cross the horizontal line at height a 
only if the above inequality holds. This directly implies that the (/-interval within 
which bound states can exist is 

r 2 _ ° ° 



-5— < g < g ( 43 ) 

q M q z — a 1 

where it is assumed that the left or the right side of these inequalities will be 
replaced by zero or infinity, respectively, when the corresponding factor r 2 — a 2 
or q 2 — a 2 vanishes or becomes negative. It also follows from that a suitable 
range of values for g always exists, so the existence of bound states in the green or 
blue region of the phase diagram is always guaranteed. Of course, the inequality 
holds also for a — > a — n -i.e., for the full right hand side of (I) (or of (II))- 
in which case the number of bound states will depend on the value of g and will 
generally be smaller than predicted by the simple positivity of the right hand sides 
of (I) or (II). Combining the above reasoning with the fact that it will always be 
q — r ~ 1 < \q\ + |r|, leads us to conclude also that e = is not a solution of the 
eigenvalue equation (I), although it can be a solution of (II) in the middle blue 
region; but ultimately even that is rejected for reasons explained in Appendix A. 
Hence the only vanishing eigenvalue is the one located in the red region and it 
represents the absolute ground state of the system. In the special case whereby the 
potential's asymptotic values Uq and Uoo are equal -i.e., when q 2 — r 2 /g- then 
clearly / m i n = f(e)\ e=q 2 — and (I) is satisfied for all positive g, a fact that also 
follows directly from the inequality (|43p . In that case, the eigenvalues are found 
easily in closed form and their formula (in the green region, say) becomes 




Figure 2: 



Phase diagram for bound states 
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in the full hypergeometric case. 

• Red region: (q > 0, r < 0). The eigenvalue e = exists and represents the absolute 
ground state of the system. 

• Green region: (q — r > 1). Bound states exist and are determined by the condition 

vV - e + \Ar 2 -ge = a- n = q- r- N (N = n + 1) (I) 

• Blue region: (q — r < for qr > or q — r < — 1 for qr < 0). Bound states exist and are 
determined by the condition 

y/q 2 — e + y/r 2 — = —a — 1 — n = r — g — n (II) 

which can be also deduced from (I) by the substitution q — 5> — q, r — > — r and iV — ► n. 
Note, however, that in the middle blue region (qr < 0) the value n = is rejected since 
then the value e = 0, which is now a solution of (II), does not correspond to a physically 
acceptable solution (see Appendix A). 

• White regions: There are no bound states. 



e " = ' 2 ~ (fevs? = " 2 ~ irf^fF"' » = » + i = m.-s ? -, (44) 

This will also hold in the blue region upon substituting q — > —q, r — > —r; while the 
further substitution N — > n is also necessary in the outer blue regions. 

Let us also mention that the wavefunction of the absolute ground state of the 
system - we say absolute since there now exist various ground states across the 
different regions of the phase diagram -is given by 

*>(*) = l^Ul • (45) 

Other features of the phase diagram worth noting are the following: 
a) In contrast to all known examples of solvable potentials, in this case the region 
of parametric space where bound states exist is doubly connected. It consists of the 
two main regions q — r > 1 and q — r < with a gap in-between (white regions) 
where no bound states exist, b) Every point (q,r), for example in the green-red 
region, has a mirror image (—q, —r) in the middle blue region with the same bound 
states except the zero one (e = 0) which exists in the green but not the blue region. 
This follows directly from the eigenvalue equations (I) and (II) and the fact that in 
the middle blue region the value n = (=> e = 0) is rejected. Thus we have pairs 
of potentials -U+ = U(q, r, x)| 9 _ r> i and U- = U(—q, —r, x)- with the same bound 
states except for e = which exists for the former but not the latter potential. 
Whether there is some symmetry behind this "pairing" of potentials in parametric 
space is an interesting question that deserves further study. 

V. CONFLUENT CASE OF THE FIRST KIND 

Let us now study the case where the parameter u vanishes, so the solutions (see 
Appendix A) are hypergeometric functions of the confluent kind -i.e., 2-P2, 1-F2 etc. 
- and their behavior at infinity will vary accordingly, including now not only powers 
of x but also exponentials. We will start with the confluent case of the first kind, 
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for which ui = but /3 ^ 0. The expression of the potential in this case emerges 
readily from ([3"3"]) and (|3~4l upon taking the limit u> — > 0, which yields 



with 





U(x) = 


Ax* + 


Bx 3 + Cx 2 4 


Vx + £ 






{x + p) 3 




A = 










B = 










C = 






+ ZpqP - q - 


^\pP + 


V = 


p 2 




-9 + \pP- 


I) 


£ = 


p 3 q 2 









(46) 



(47) 



and has a typical plot as in Figure 3. 

t U(z) 




Figure 3: Typical potential of the confluent family of the first kind. At — oo the potential 
tends to a fixed value and at +oo it rises as a harmonic oscillator. 

Let us note here that some of the results for the confluent case derive directly from 
the previous ones by setting oj — 0, provided they do not represent a qualitative 
change in the problem, in which case it is safer to redo the calculation. Let us 
simply mention these results without elaborating much on them. We begin from 
the relation z — z(x), for which we have 



z , {x) = V 1 + Wp) ^ x{z) _^ e , ; x{z) _^ p i (48) 

X — z— H-oo 4 

The general expression of the wavefunctions is now 



iK*) = X (x 2 , ^yt 2 ( X ( X + P) F ' + ( aX + P( a + ^) F ) > ( 49 ) 



(x + p) 

where F(x) = 1F1 (a, 6; c, d] —j3x) and 



a = -q + \Jq 2 - e c = 2 - q + ^ q 2 -e (= a + 2) 



b = -q+ ^Jq 2 -e + pe+1 d = 1 + 2yV 



(50) 
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with q as before -formula (30 a )- and 

where the new parameter p takes now the position of r in the phase diagram that 
we will present shortly. The search for physically acceptable solutions depends 
crucially on the sign of /3, or p, and is based on the asymptotic relation 

k->oo I (6)1 (c — a)I (a — a) I (a)l (c — 6)1 (a — o) 

T{a)T(b) 

The results we obtain are as follows. 

A. Scattering states 

Since the only scattering that makes sense now is from the left, we need only 
calculate the reflection amplitude r^ik) = r(k), which becomes 

m ~ q + ik T(l-2ik)T(l-q + pe + ik) { P) {l < ° J 

(53) 

r(k) - -1 + lk r(l + 2ik)T(q-pe-ik) 2lk 

r{k) - T+ik T{l-2zk)T{q~pe + ikf {P > °> 

with \r(k)\ = 1 of course, since the form of the potential does not allow for the 
particle to pass through towards positive infinity. 

B. Bound states 

A similar calculation as before yields now -apart from the vanishing eigenvalue 
in the region q > 0, p < 0, which is found separately (see Appendix A)- the two 
complementary conditions 



y/q 2 - e = -pe + q - N (p<0,N^l) (54) 

and 

\/ q 2 — e = pe — q — n (p > 0, n ^ 0) (55) 

from which we can readily obtain the results summarized in the next phase diagram. 

Clearly, the phase diagram has the type of mirror symmetry we mentioned before 
and this raises similar questions as to its origin and interpretation. 

Concerning the absolute ground state of the system -i.e. the one with e = 0- 
this exists only in the red region and its wavefunction is written as 

x q e /3x/2 

M') = j^pi- (56) 

Note also that the bound states do not depend on the full set of dimensionless 
parameters p, /3, q of the potential U(x) (in the invariant parametrization as before) 
but only on their two combinations p(= l/p/3) and q. This is contrary to what 
happens in the scattering data -formulas (|53[) - where all three parameters enter. 
Something similar to this occurs in the full hypergeometric case we examined before. 
The presence of isospectral orbits in parametric space is thus one more feature of 
the potentials presented here which has to be explained. 
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Figure 4: Phase diagram for bound states 
in the confluent case of the first kind. 

• Red region: (q > 0,p < 0). The eigenvalue e = exists and represents the absolute 
ground state of the system. 

• Green region: (p < : p ^ (q — l)/q ). Bound states exist, and they are given by the 
formula 



£n = 



{pq-h) ~PN + sj \pq-\) +pN 

p2 



N = 1,2, 



1= (i -pq)q 



(I) 



Blue region: (p > : p l/q for pq > and p ^ (q + l)/q 2 for pq < 0). Bound states 
exist, and they are given by the formula 



(pq - |) +pn + \/ [pq - |) 2 —pn 



p 



n = 0, 1, 



(P9 - 1)3 



(II) 



which can be also derived from (I) by the substitution p — > —p, q — > — q and iV — > n. 
But in the left blue region (pq < 0) the value n = is rejected for the same reasons as 
before. 

White regions: No bound states exist. 



VI. CONFLUENT CASE OF THE SECOND KIND 

The expression for the potential U(x), in the invariant parametrization, emerges 
directly from (|T7|) upon setting /3 = 0, so there will be 

m = ± ^ — A-y, — 1 ( 57 ) 

with x = x(z) as in (|48|). which implies an exponentially fast approach to the 
asymptotic value Uq = q 2 of the potential at z — > — oo, but a very slow decrease to 
zero -of the kind U ~ 1/z 2 - at z — > +oo. Depending on the range of values of q 
the potential U(z) can take one of the following forms: 
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Figure 5: Representative forms of the potential U(z) 
in the confluent case of the second kind. 



The wave functions will now be written as (Eq. (|49| with /3 = 0) 



1>(x) = y/4 « X + P) F ' + ( ax + P( a + V) F ) ( 58 ) 

where F(x) — iF2(a;b,c;—ex/p) and 

a = -q + vV - e, 6 = a + 2, c = 1 + 2vV - « (59) 

while the asymptotic behavior of -F(£) for large |£| -where £ = — ex/ p in our case- 
will be given by 

F(n > r(b)r(c) 

^ ~KhW r(6-a)r(c-a) 1 U 

+ M«Ke<<™«^ + e-<<*^) (■» " K- c + i)) (6o) 

e.g., by a linear combination of the three asymptotic behaviors at infinity, which 
now have the form of two exponentials and one power of £ due to the fact that the 
maximum dimension component of the operator C in (f2"2")l -with u> = f3 = 0- is 
now of first order. 

On the basis of the above the results we obtain are as follows. 

For e > Uq there are scattering solutions both from the left and from the right, 
and the corresponding reflection amplitudes are given by the formulas 

r,(k)= -^i k If + ^J -Y™*-"*, r R (k) = ie-^ (61) 
q + ik I (1 — 2iK) \p J 

where it is noteworthy that, apart from a constant phase (independent of k), the 
reflection amplitude from the right is real. 

As for bound states, there is only the eigenvalue e = with a corresponding 
eigenfunction 
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which will satisfy the boundary conditions ip(0) — ip(oo) = only if < q < 1/4. 
Nevertheless, (|62p is not square integrable as one might have expected, due to the 
position of the corresponding eigenvalue at the threshold of the continous spectrum. 

It hardly needs mentioning, that all the above results -irrespective of how 
they were produced- can be verified directly with a straight substitution in the 
Schrodinger equation. To do that, it is advisable to cast that equation in the equiv- 
alent x-form 

where the derivatives and the expressions of the functions w and U are with respect 
to the initial variable x. 

Finally, we note that in both main classes of potentials (full hypergeometric case 
and confluent of the first kind) their graphs do not always assume the typical forms 
of Figures 1 and 3 but they show interesting variations similar to those in Figure 
5. The question of which of these forms may have some special physical interest 
will not be discussed here. 

VII. PENDING ISSUES AND GENERALIZATIONS 

Let us note first that our preceding discussion on potentials that are solvable via 
hypergeometric functions of third order is not exhaustive. The first reason for this 
is that the solution yo — x s (x + p) , which we have chosen to reduce the order of the 
initial equation, is just a special case -for s' = 1- of the general expression, tjq — 
x s (x + p) s , for the solutions of the first order equation My — 0. This special case 
is undoubtedly the most important, for reasons that may have become obvious by 
now. But there are also other possibilities -e.g., s' = 2- belonging to the terminating 
series category and which can thus be solutions of Ly = 0, as required. Given 
though that the coefficients of the relevant polynomial are now related to one 
another (they are all functions of p) , the number of independent parameters of the 
operator L, and hence of the potential V, will necessarily be reduced, or such a 
solution might not even exist. 

But there is another reason why our previous discussion is incomplete. We as- 
sumed that p and ui are positive -and hence, they fall outside the region < x < oo 
between the singular points of the original equation- which means that the poten- 
tial V is also finite in this region, and therefore also in — oo < z < +oo into which 
< x < oo is mapped via the transformation z = z(x). But if we assume, for 
example, that uj < -and substitute u) for —ui— then x must lie within the inter- 
val < x < which is mapped via the transformation (|23j) to the seminfinite 
z-interval — oo < z < 0; and so the corresponding potential V(z) will also be a 
half-interval potential with repulsive core of the type 1/z 2 , for z — > 0. In the same 
fashion one needs to examine all the remaining choices of signs for p and uj together 
with the corresponding choices of intervals between singular points with respect to 
x or z. What is certain, is that even in the context of the third-order case, the 
set of solvable potentials is much wider than what we presented here and deserves 
further study, especially if such a study could provide answers to some interesting 
special questions. 

A third pending issue pertains to the comparison of the present family of solvable 
potentials with that of Natanzon. What are the differences or similarities between 
the two families? Can one family be deduced from the other by applying the special 



17 



Darboux transformation related to the removal of the ground state? We will leave 
this discussion for a follow-up publication, once we have constructed a much more 
general framework of studying solvable potentials and developed a general formal- 
ism that is immediately applicable to Natanzon potentials as the simplest special 
case. What is certain is that -apart from the special case p = lo^ 1 (for which the 
potentials are elementary functions)- there is no overlap between the two families, 
as evidenced both from the different scattering data and the different functional 
forms of their solutions. (Hypergeometric functions of the type 3F2, 2-F2, 1-F2 m our 
case, vs. 2-F1, in the case of Natanzon.) 

Next comes the question whether our methodology can be extended to equations 
of order higher than third, and produce the corresponding families of solvable 
potentials. While this will be the topic of an upcoming publication, we find it 
warranted at this point to present, without proof, some of the pertinent conclusions 
in order to underscore the existence of a general framework for the systematic 
investigation and enlisting of, potentially, all solvable potentials. Our most basic 
results are contained in the following three propositions: 

Proposition 1: Every formal eigenvalue equation of the general form 

(L + \M)y = (64) 

where L is any linear differential operator of order n and M a similar operator of or- 
der 

n — 2, can be reduced to an eigenvalue equation of second order -and hence, to 
a Schrddinger equation- provided that each solution of My — is also a solution 
ofLy = 0. 

Proposition 2: The weight function w{x) of the Liouville equation that results 
from the reduction of any given bidimensional equation (L + XM)y = -subject to 
the constraints of Proposition 1- can always be written in the form 

X z (l + LOX) 

Therefore, the transformation function z(x) that renders the Liouville equation into 
Schrddinger form will also be a "universal function", i.e., the same function for all 
solvable potentials. 

Proposition 3: The solvable potentials arising from a bidimensional equation of 
an arbitrary order (L + \M)y = -subject to the constraints of Proposition 1- will 
always have the functional form 

T _, , Ax i + Bx 3 + Cx 2 +Dx + E 

v ( x ) = 7 ^rr, n (66) 

v ; {x + pf{l + wx) v ' 

where the coefficients A, B, C etc., are functions of the parameters that remain 
free in the operators L and M once we demand that every solution of My = be 
also a solution of Ly = 0. 

In subsequent publications, we will rely on these three propositions -and on some 
crucial extensions of them- in order to extend this work and develop a systematic 
theory of exact solvability of the Schrodinger equation. 



VIII. DISCUSSION 
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As was noted at the outset, the critical new idea of this work is this: Given 
that the Natanzon search for exactly solvable potentials is exhaustive, we can only 
produce new ones if we enlarge the set of functions within which their solutions are 
sought for. And if this extended set were the hyper geometric functions of all orders, 
then utilizing this set for solving Schrodinger equation would be in principle possible 
only in conjunction to the idea of order-reduction. We are thus led naturally to the 
notion of eigenvalue equations of higher order -of the hypergeometric type if we 
want them to be solvable- that can be reduced to the Schrodinger equation using 
the technique of order-reduction. And as we saw earlier, this idea is indeed feasible 
and has actually delivered as a tangible product a new set of solvable potentials, 
beyond that of Natanzon. We also saw that this notion -provided that propositions 
1 through 3 are valid- is extensible to a very elegant general formulation which 
encompasses, as special cases, both the one studied here -where L was of third and 
M of first order- as well as the case of Natanzon, with L being of second and M 
of zeroth order, i.e., a function of x. 

It follows that the method is in principle extensible to equations of higher than 
third order, even though it remains to be seen whether the relevant calculations 
continue to be manageable. Many further questions arise; here is but a short list of 
these: 

a) Could there exist an infinite hierarchy -an infinite "tower"- of solvable po- 
tentials, with those of the Natanzon class occupying its "ground level", the ones we 
presented here the "first floor" and so on, for the potentials originating from equa- 
tions of ever higher order? Or does the process terminate leading to an ultimate 
closed set of solvable potentials? 

b) What is the origin of the peculiar topology of the phase diagram for bound 
states and its mirror symmetry! Are the potential pairs (7+ and U— , connected by 
this symmetry, supersymmetric partners? And what about the property of shape 
invariance that was shown not to match the initial expectations pertaining to its 
range of applicability? In view of the preceding discussion, could it be that an old 
conjecture 3 , namely that this invariance may require a wider class of potentials 
than that of Natanzon, ought to be investigated further? 

c) What is the origin of the reparametrization that reduces the number of active 
parameters of the system, or of the isospectral orbits in parametric space? Do 
these properties hint at some sort of symmetry that allows also a purely algebraic 
approach to the problem? 

At a more mathematical level, the mapping of eigenvalue equations of higher 
order into the Schrodinger eigenvalue equation -a mapping that can be cast in an 
elegant compact form- could potentially provide valuable new insights to an old 
equation which clearly continues to hold an element of surprise. 

APPENDIX A 

CALCULATION OF EIGENVALUES AND EIGENFUNCTIONS 

According to the equations ([T]) through (JTUJ) of section II, we have 



ip(x) = (x + p) 7 ^x q (l + Lux) r - q+ ^^ 



( 



x s (x + p) 



II 



) 



(Al) 
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where q and r are defined as in Q30p and y(x) is the solution of the bidimensional 
equation (|22p . The latter can be written in terms of the hypergeometric function 
3F2 as 

y{x) = 3 F 2 (a,b,c;d,e;-uJx) (A2) 

where the factor —ui, in the argument £ = —cox, serves to transfer the finite singular 
point xq = —lu^ 1 of (|2"2"j) to the standard position (b = 1 of the hypergeometric 
equation, while the factor x^ represents the physically acceptable power behavior of 
© for 

x —> 0. The possible values of /i (physically acceptable or not) are found as so- 
lutions of the unidimensional equation (or Euler-type equation) 

Cix" = (A3) 

where C\ is the lowest-dimension component of the bidimensional operator C = 
L + AM of (|2"2")l . For the power behavior x v at infinity the respective equation is 

C 2 x u = (A4) 

where C 2 is the largest-dimension component of the operator C. (A3) and (A4) are 
of course cubic equations for /1 and v, which are however easily solved since we 
already know (from the known solution yo = x s (x + p)) that one behavior at the 
origin is x s and one behavior at infinity is x s+1 . The exponents of the remaining 
two behaviors are easily found to be 

— 3 — a — s , r— _ 2 — (fi/u)) — s 1 /— 

2 2 

where the specific choice of signs for /1 and v (positive for /1 and negative for v) 
reflects the fact that the physically desirable behaviors are then a; M and x v for 
small and large x respectively. Now we can readily calculate the hypergeometric 
parameters a,b,c,d and e in (A2), by recalling that the parameters a$ and bj in 
the general hypergeometric series 

C / t 1 1 \ (<Zl)ra • ■ • (Op)n Z n / * fj\ 

p F g (oi,a 2 , • • .,a p ;b 1 ,b 2 ,.-.,b q ;z) = 2^ r (A6) 

have always the same asymptotic meaning: x~ ai (i = 1, . . .,p) are the power be- 
haviors of the general solution of the corresponding hypergeometric equation at 
infinity and x°, x 1 ^ bj (j = 1, . . . ,q) are the power behaviors at the origin. The 
calculation of a, b, c, d and e in (A2) reduces then to simple asymptotic comparison 
of the two sides for x — > and x — > 00. The result is 

a = u — s — 1 , b = u, — v. c = u — V 

(A7) 

d = fi — s + 1 , e = /i — + 1 

from which, upon using (A5) and the definition of the new parameters q and r (eq. 
(f3T))0 . we obtain the expressions (|3l)t and the general formula (f3l>f for the solution 
i/>(ir) that satisfies the boundary condition ^>(0) = in the region of bound states. 

The boundary condition at infinity (?/>(oo) = 0) is imposed via the asymptotic 
form of ip( x ) f° r large x which is written as 



ip(x) -> (a-b)Bx-^-^ +{a~c)Cx^~^ (A8) 
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and is derived from the asymptotic relation (f3T)|) for F(x) 

F(x) = 3 F 2 -> Ax' a + Bx- b + Cx~ c (A9) 

by taking into account what we noted earlier: Namely that the first-order operator 
M. acting on F (Eq. ([35])) annihilates the term x~ a . It follows then from (A8) that 
the boundary condition at infinity will only be satisfied if a — c = or C = 0. But 
the case a — c = is rejected at once due to the factor T(a — c) in the expression 
for C, so we are left with C = which can be realized in each of the following ways 

i) a — —n ii) b = — n iii) e — c = —n (A10) 

since the case d — c = — n is also rejected due to the relation d = a + 2 and the 
presense of the factor T(a — c) in the numerator of C. From T(d) /T(a) — a(a + 1) it 
is also clear that for case (i) only the values a = and a = — 1 survive, the former 
reproducing the already known ground-state solution e = (i.e., A = Ao) for q > 0, 
r < 0, while the latter gives i/j = as expected. In case (ii) the pertinent condition 
is written as 

Vq 2 - e + Vr 2 - ge = q - r - (A =1,2,...) (All) 

and will be satisfied under the conditions laid out in section IV (Eq. ((H])) from 
where the exact number of bound states for the problem is deduced. Note that 
e = is no solution of (All) for n = 0, i.e., for N = 1, and thus the condition 
(ii) yields only positive eigenvalues. In both the above cases, the hypergeometric 
function F(x) has polynomial form because both conditions a = — n and b = —n 
are termination conditions for the respective hypergeometric series. 

However, in case (iii) -contrary to what usually happens for bound states- the 
hypergeometric series does not terminate since in that case both behaviors x~ a and 
x~ b survive at infinity, even though only the latter contributes to the asymptotic 
limit of ip(x). The condition for calculating the eigenvalues is now written as 

V<1 2 - e + Vr 2 - ge = r - q - n (n = 0, 1, 2, . . . , ^ r - q) (A12) 

but especially in the middle blue region of the relevant phase diagram (Fig. 2) 
the value n — has to be exempted since then (A12) can have as its solution the 
vanishing eigenvalue e = 0, which, however, does not lead to a physically acceptable 
solution tp(x) and must be thus rejected. The reason for this is nontrivial. For e = 
a flaw arises in the mechanism that "ejects" the undesirable asymptotic behavior 
for x —> oo, since then the hypergeometric parameters a, b, etc., no longer depend 
on p (see Eq. ([351)1. so the corresponding solution F(x) has no way of "knowing" 
the value of p although the relevant ejection operator Ai — x(x + p)d + ax + p(a + 1) 
in (|35p continues to depend on it. 

Using the same techniques as above -but with the asymptotic formula ([5^]) in- 
stead of ([Sn])- we can also calculate the bound states in the confluent case, obtaining 
as a result -except for the zero eigenvalue- the two conditions 

\/ q 2 - e = — pi. + q - N, \/ q 2 — e = pe - q — n (A13) 
which lead to formulas (I) and (II) of Figure 4. 
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As for the scattering solutions (in the full hypergeometric case for example), we 
need the asymptotic form of ip(z) for z — > ±00 which turns out to be 

2-^-00 pi/i[a + l)e lkz 
^(z) = / (A14) 
\ z— »oo ~Be~ ik ' z + ~Ce ik ' z 

(B = uj—^/^ia - b)B, C = co—^ia - c)C) 

and it is clear from (A14) that the scattering solution from the right equals ifj* 
-that is, ipR = ip*- and therefore 

(A15) 

from where -after a few steps of algebra- we readily obtain the results (|4"U| and 
(|4"T1) . Finally, the calculation of scattering from the left is done by forming the 
combination c\i\) + C21P* containing only a traveling wave exp(ifc'z) at +00, and 
proceeding as before. The spectral data for the confluent cases of first and second 
kind are calculated in an analogous manner. 

Some closing remarks about the connection of hypergeometric functions and 
bidimensional equations are in order. The connection originates from the fact that 
bidimensional equations lead to two-term recursion relations, from which the gen- 
eral coefficient of the series-solution is computed in closed form and it turns out 
to be exactly the same as the general coefficient of the hypergeometric series (A6) . 
Thus the hypergeometric functions emerge naturally as solutions of the bidimen- 
sional equations of step unity. The equation they satisfy -take, for instance, the 
function 3F2- will naturally have the form ([T5]) . but with oj = — 1 (so that the 
nonzero singular point lies in the standard position xq — 1) and e = 0, so that 
one exponent of the power behavior at the origin can take the typical value zero 
as in the general hypergeometric series (A6) that solves the relevant equation. As 
for the exact form of the remaining numerical coefficients of (TT5|) . as functions of 
the hypergeometric parameters a, &, c, d and e, this is never needed since to obtain 
the solutions it suffices to know the power behaviors of y(x) at zero and infinity, 
together with the standard asymptotic meaning of the hypergeometric parameters. 

As for the hypergeometric functions of the confluent type -i.e, with p < q + 1- 
these emerge naturally from those bidimensional equations whose highest-dimension 
component L2 has no longer the full order of the equation, but is of lower order by 
one, two, three, etc., exactly as the value of p compared to q+1. In this case, only 
p behaviors at infinity -where p now stands also for the order of L2- will be powers 
of x; the rest will be exponentials. This alters the character of the point at infinity 
that now becomes irregular singular. As for the terminology, we have chosen the 
terms confluent equation (or function) of the first kind if p = q, confluent of the 
second kind if p = q— 1, etc. To complete the picture, we note that not only the 
bidimensional equations of step one can be solved via hypergeometric functions, but 
also the bidimensional equations of any step £. The reason for this is simple. It has 
to do with the (rather obvious) fact that the transformations t = x m (change of 
independent variable) and y = x^Y (change of dependent variable) preserve the 
bidimensional character of an equation and merely cause a change in its step, from 
£ to £/m, and a shift in its starting powers by /i, respectively. 
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It thus follows that we will always be able to write the solutions of any bidimen- 
sional equation in the form 

y(x) = x^ p $q(ai, . . . , a p ;bi, . . . ,b q ; kx e ) (A16) 

where p $ q is the general solution of the hypergeometric equation with parameters 
aj and bj, £ is the step of the given equation and k the numerical coefficient that 
maps all singular points outside the origin (all of which lie on a circle in the complex 
plane) to the standard position ( = 1 where £ = kx l . As for the parameters \x and 
di, bj in (A16) these are obtained readily from the asymptotic comparison of its 
sides at zero and at infinity, leading directly to the relations 

Vi = n-£ai, Hj = fj, + £(1 - bj) (A17) 

where jij and vi are the power behaviors of the given equation at zero and infinity 
respectively, while of course \x must be equal to one of the [ij (e.g. = fi\). 

It is clear from the above discussion that the hypergeometric functions and 
bidimensional equations emerge as the natural framework for the systematic study 
of the problem of exact solvability of the Schrodinger equation or any other equation 
for that matter. Our future work on this problem will lie within this systematic 
framework. 
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